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1. INTRODUCTION 

Low density parity check (LDPC) codes [1] have become increasingly popular, due to their 
excellent performance and efficiency, to the extent that they are being considered as the error correcting 
coding technique of choice for major standard developing bodies such as the third-generation partnership 
project (3GPP) and the IEEE. Among various decoding methods, decoders based on belief propagation [2], 
[3], offer superior performance comparing to other hard decision decoding algorithms such as bit flipping 
(BF). However, their optimal performance is very computationally intensive. On the other hand, decoding 
techniques such as BF method [1] are easily implemented and offer relatively low computational complexity. 

BF algorithms have been extensively studied and coding researchers have offered many variants of 
the method including weighted bit flipping (WBF) method [4] where soft channel information is 
incorporated, and a reliability measure is introduced to enhance performance of the BF algorithm. More 
recently, Wadayama et al. [5] devised a class of BF decoding method, i.e., gradient descent bit flipping 
(GDBF) algorithm that is based on maximum likelihood (ML) decoding problem. The GDBF method follows 
a gradient descent optimization model for the ML decoding to maximize an objective function [5]. 
Sundararajan et al. [6] introduced a noise perturbation to the objective function in order to avoid being 
trapped in local extreme, that do not correspond to a solution, during decoding process. 
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GDBF and noisy gradient descent bit flipping (NGDBF) algorithms are considered probabilistic bit 
flipping methods, and from that point of view, a few other methods have been proposed and implemented to 
improve decoding processes. Rasheed et al. [7], for instance, introduced a random perturbation to the 
decoding scheme, and subsequently, a random set of bits (that minimizes the inversion function) are selected 
to be flipped. The probabilistic gradient descent bit flipping (PGDBF) method proposed in [7], however, is 
designed for binary symmetric channels (BSC). A few other improvements have also been made, such as the 
work of Cui and Wang [8], in which a tabu-list is designed to overcome the repeated (the same) bit flipping 
issue. In other words, the set of flipped bits in the current iteration are added to the tabu list to prevent them 
from being flipped in the next iteration. Cui and Wang [8] improved their method in [9] by proposing an 
information storage bit flipping (ISBF) decoder for BSC channels in which the sum of syndromes orthogonal 
to a bit (also known as energy of the bit) is calculated. Then a bit is flipped with a certain probability if its 
energy value reaches a certain threshold (Th). By storing the energy value from the previous iteration, the Th 
can be set as the maximum of both energy values. In this method, both threshold and energy values can be 
computed in parallel, therefore reducing decoding delays. 

The class of NGDBF method initially proposed in [6] as well as other variants of the method so far 
provide significant improvements over GDBF method for transmissions over additive white Gaussian noise 
(AWGN) or BSC channel. In fact, Declercq et al. [10] have shown that introducing randomness (noise) to an 
otherwise noiseless decoder can greatly improve error performance of family of BF decoders (including 
NGDBF and PGDBF algorithms) such that their performance can be compared with soft decision decoders’ 
such as sum product algorithm (SPA) or other variants of belief propagation decoders. 

While performance of the NGDBF algorithm and its variants show vast improvements over regular 
GDBF algorithm there are two major issues associated with them; to achieve performances similar to soft 
decision algorithms, NGDBF method require large number of decoding iterations [10] which may impose 
practical limitations on applying the decoders to use cases where low latency is required. It is because end to 
end latency of the system is directly related to processing delay associated with decoding iterations. Second, 
NGDBF algorithm, as it is currently designed, performs well over AWGN channel. However, its 
performance may degrade where it is to function over more realistic channels such as Rayleigh fading 
channel over which a lot of more modern use cases and applications are envisioned to operate. 

In this paper, an improved noisy GDBF method is proposed that efficiently decodes LDPC encoded 
codewords transmitted over correlated Rayleigh fading channel. We demonstrate in section 3 that our 
proposed method outperforms regular NGDBF decoder in 1) average number of iterations required to achieve 
a certain probability of bit error, and 2) error performance in the form of probability of bit error as a function 


of energy of the bit per noise power spectral density (fo J. No" In section 3, we explain our proposed method 


chronologically and follow through each stage of a communication system as bits are encoded, transmitted 
through the channel, and decoded using our proposed and improved noisy gradient descent bit flipping 
algorithm. The results of our simulations are presented in section 4. 


2. RESEARCH METHOD 
2.1. Encoding and channel model 

Assuming binary phase shift keying (BPSK) transmission, a codeword V = (Vo, V1, V2, .-,UVn—1) iS 
modulated to a sequence x = (Xo, X1, X2, -.-, Xn-1). For our simulation, the generator matrix g(x) = 1+ x? + 
x*4+x°+x7+4x'° is used to generate codeword v resulting a (21,11) projective geometry low density 
parity check (F) code. The modulation scheme maps each bit x, to vz where x, = 2v, — 1 with k = 
0,1,...,2 — 1. If x is transmitted through an AWGN channel the received vector Yawgn is of the form 
Yawgn = (Yidwgris i erro where Yrawgn = Xx +N, with ny, being a random variable that is 
Gaussian and has zero mean and variance Ny/2 (noise power spectral density). Assuming a binary operation, 
Xg = 2Vk — 1 = +1 for vy, = 1 and x, = —1 for v, = 0. Subsequently, soft decision received bit Yk awgn = 
Xg +n, =1+n, and —1 + ng, respectively. On the other hand, if x is transmitted through a Rayleigh (flat) 
fading channel, the received vector Yyay is of the form Y;ay = (Visas ry Vi eeay ys where Yk ray = RexXK + 
Ng, with Rg being a zero-mean complex Gaussian random variable. The absolute value of Rẹ (the envelope 
of the channel response) has Rayleigh probability distribution [11]. And finally, z = (Zp, Z4» Z2,-.-,Zn—1) is 
the hard decision received vector estimated from y vector. In our simulations, we let z; = 1 if y; = 0 and 
Zi = 0 otherwise. 

It has been shown that discrete time samples of a Rayleigh fading process must necessarily be 
correlated and the correlation function is dependent on the Doppler frequency associated with relative motion 
of the transmitter, receiver, and other parameters such as propagation path and antenna characteristics [12]. If 
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a very long bit interleaver is employed (where bits are spread over an interval much longer than the channel 
coherence time) then R,'s are independent (and complex) random variables [13]. However, in case of no (or 
partial) interleaving, R;,’s become more correlated random variables which may introduce blocks of bit error 
in the transmitted frames when the transmitted frame experiences a deep fade. 

We used the young model to simulate the channel coefficients R;,’s that is based on the inverse 
discrete Fourier transform (IDFT) technique described in [12]. We generated time samples of Rayleigh 
fading processes corresponding to different Doppler frequencies (shifts) to simulate Rayleigh fading channel. 
To ensure that appropriate channel coefficients are generated, we calculated the mean of the generated 
random variables R, and observed it to be zero or approximately zero. It was in addition to verifying level 
crossing rates that vary as a function of Doppler frequency in a Rayleigh fading channel model. To quantify 
the extent of correlation in a given Rayleigh channel model, we chose a comparative method in which the 
channel coefficients are compared to the receiver noise that is independent and identically distributed, i.e., 
the normalized cross correlation of Rg’s and n,’s in the received signal Yg ray- 


2.2. Decoding 

BF decoding, in its simplest form, is a syndrome decoding technique where parity check matrix of a 
code H determines the syndrome of a code s in error detection and correction as described below. Assuming 
H to have J rows and n columns, the number of 1’s in each row and each column of the parity check matrix 
are the row and column weight of the matrix and are denoted p andy respectively. For an LDPC code, 
pandy are small comparing to J and n, hence the name low density. The syndromes [1] of the received 
vector s is then defined: 


S = (S1, Sz, =, SJ) = Z. HT (1) 


where z is the hard decision received vector described in sub section 2.1 and the superscript T denotes 
transpose of the matrix H. When errors are introduced in the received vector, there will be parity failures and 
some of the syndrome components become 1. BF algorithm calculates all parity check equations (check 
sums) and changes (flips) the bit(s) that are contained in greater than a set number 0 of unsatisfied parity 
check equations. It has been shown that @ (also referred to as the threshold) is a function of the p,y, error 
correcting capability of the code and the signal-to-noise ratio (SNR) [14]. 

BF algorithm first calculates the syndrome bits (parity check equations). If all syndrome bits are 
zero, then the process stops and the correct codeword is declared. Otherwise, it finds the number of non-zero 
check sum equations for each bit k, denoted f,, where k = 0,1, ...,n — 1. The BF algorithm then identifies 
the set of bits for which fẹ is maximum and flips the bits in the set. The process is then repeated until all 
check sums are zero. To reduce processing delay, this iterative process may be limited, that is the process 
stops if all check sums are zero or the maximum number of iterations is reached. 

To improve hard decision BF decoding, some measure of reliability is introduced to the decoding 
process. For an AWGN channel, the magnitude of the soft decision received bit |y,| may be used. The larger 
|¥;| the more trustworthy the hard decision bit z, will be. Kolesnik [4] used this concept and developed an 
improved BF decoding method known as weighted BF (WBF) decoding. In WBF algorithm, a new function 
E,, is defined for each bit k: 


k 
E, = Ewer — 11/1 n (2) 


where Sẹ is the set of check sum equations that contain bit k (bit k participates in them). And | Meee is the 


magnitude of the least reliable soft decision bit. A bit is flipped whose Ep (or energy function) is greater than 
zero [4]. A few modifications and improvements have been made to WBF algorithm. WBF algorithm 
determines whether a bit should be flipped based on the contributions from the check sums as shown in (2). 
Zhang and Fossorier [15] proposed a modification to (2) by incorporating a factor such that the reliability of 
the individual bits also help determine if a bit should be flipped: 


1 k 
E', = Ewer -DIYO in Z Yk (3) 


where æ is areal number and a > 0. The value of parameter æ was found to be a function of SNR and the 
column weight y of the code, regardless of the type of the LDPC code. As for the flipping threshold of the 
family of BF decoders, Cho and Sung [16] proposed a method that adaptively adjusted the flipping threshold 
in order to avoid meaningless no flipping conditions and to improve rate of decoding convergence to an 
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optimal solution. The flipping threshold was shown in [16] to be a function of the number of errors. To avoid 
no flipping condition, the maximum of all E’, in (3) Emax is set as the new threshold. Initially, any bits 
whose energy value is larger than the initial threshold is flipped. In case of no flipping condition, the 
threshold is adjusted to Emax and the process is repeated. 

Wadayama et al. [5] proposed the original GDBF algorithm in which the maximum likelihood 
decoding concept presents itself as an objective function for gradient descent optimization. Based on 
maximum likelihood (ML) decoding principle, the received bit that maximizes the objective function f(x), 
defined in (4), is declared as the correctly transmitted code bit. 


F%) = DRO XK Ye — Ly te5, (25/ — 1) (4) 
J 


The first term in (4) is the sum of the cross correlation of the transmitted code bits and the soft decision 
received bits which needs to be maximized to correctly decode the transmitted codeword. The second term 
plays the role of a penalty term. If all check sums are satisfied, the second term has an absolute minimum of 
—p where p is the row weight of the parity check matrix. To ensure the objective function achieves its global 
maximum at the end of the iterative decoding process, partial derivative of objective function F,, defined in 
(5) is evaluated bit by bit via a gradient descent optimization process. 


Ex = Xk: Yk — Des -1) (5) 


As for mechanics of bit flipping, there are couple of points to consider; the first point is the fact that 
flipping bits can be deterministic or probabilistic. In the former, a bit that meets certain condition is flipped. 
In case of BF (and its variants) algorithms, the condition is number of unsatisfied parities checks that are 
orthogonal to a bit. In the latter, once the condition is met, the bit is flipped with certain probability. One 
study [17] has shown that the probabilistic BF algorithm is more useful when the number of unsatisfied 


parity check equations is close to [4] with J being the number of rows of the parity check matrix. They 


showed that for a binary symmetric channel, the probability p with which bits are flipped is a design 
parameter and depends on the cross over probability of the channel. In general, the greater the number of 
unsatisfied parity check equations the higher the flipping probability p. The second point regarding the 
mechanics of bit flipping, is the number of bits that can simultaneously be flipped. From that perspective, 
there are three flipping modes: single bit flipping, multi-bit flipping, and a combination of both (or hybrid). 
In case of GDBF decoding where an optimization process is ensued based on gradient descent (or ascent) 
method, bit flipping can be considered as steps on a search path toward the solution that is the correctly 
decoded codeword. Single bit flipping is equivalent to taking small steps on the search path that may 
converge slowly to the solution. Multi bit flipping, on the other hand, is equivalent to taking larger steps on 
the search path. The hybrid bit flipping (also referred to as mode switching) usually starts with multi bit 
flipping and switches to single bit flipping as the search approaches the solution [5]. To determine how many 
bits are to be flipped in a multi bit flipping process, a threshold @ is defined for inversion function defined in 
(5), and the bits whose inversion function fall below the threshold @ are flipped. The hybrid and multi bit 
flipping methods have been shown to have better performances comparing to single bit flipping described 
here [5], [6], [18]. 

The objective function defined in (4) is not a linear function and may have many local maxima in 
addition to the local maximum that corresponds to the transmitted codeword. During the search for the 
transmitted codeword, it is possible that the search point may be trapped on a local maximum that does not 
correspond to the transmitted codeword. One method to escape the so-called trap, proposed in [5], is to add a 
perturbation to a trapped search by switching from single bit flipping to multi bit flipping described above. A 
disadvantage to this method is that a multi bit flipping process may miss the desired solution due to the fact 
that the search steps may be too large. Another method to escape the traps, proposed in [6], is by adding a 
random perturbation to the inversion function (5): 


Eg = Xk: Yk — (w Oes — 1) + qk (6) 


: i : i : No... 
where w € Rt and qx is a Gaussian random variable with zero mean and variance a? = q? ra withO<as 


1. The study [6] found that w and q@ are close to 1, and their optimal values are code dependent, however, 
they are weakly dependent on signal to noise ratio. The decoding method proposed in [6] is referred to as 
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noisy gradient descent bit flipping (NGDBF) algorithm. Single bit flipping and multi bit flipping techniques 
can both be performed. The study also used adaptive threshold technique when multi bit flipping method was 
used. As far as error performance of NGDBF method is concerned, the method outperformed the regular 
GDBF method in almost all scenarios, some of them, however, at the expense of more iterations. A more 
recent study [19] added more solid theoretical framework to GDBF and NGDBEF algorithms that had been 
otherwise empirically analyzed. In the study, the authors proposed a decoding technique that is based on the 
conditional probability of message error given the syndrome information, defined in (6). The decoding 
technique proposed in [19] is also based on maximum likelihood decoding principle and is referred to as 
probabilistic local maximum likelihood bit flipping (PLMLBF) algorithm. 


2.3. Improved NGDBF algorithm 

To the best of our knowledge, none of the decoding algorithms described in this paper, considered 
real life applications where AWGN channel may not appropriately model the channel over which codewords 
are transmitted. Some of the more modern applications in wireless communications realistically occur where 
other channel models such as Rayleigh fading channel may be better suited to model the communication 
channel. In those cases, the error performance of BF decoder and its variants such as NGDBF may be 
degraded. Additionally, while decoding in non-Gaussian channels such as Rayleigh fading channel, BF 
method and its variants may suffer large processing delays due to the fact that they may have to decode using 
greater number of iterations than normal. To put this issue in perspective, one application is here described in 
the context of intelligent transportation system (ITS), i.e., vehicle to vehicle (V2V) communications. 

V2V communication enables vehicles to wirelessly exchange information about their speed, location 
and heading in order to avoid crashes. The information is broadcast by each vehicle and is received by all 
other vehicles in a given range. Given crash avoidance is a safety critical application, V2V communication 
must have high reliability and low latency. The required latency for V2V communication is application 
dependent and ranges from 10-100 ms [20]. Many V2V applications, described in [20] and [21], may arise 
while vehicles are moving at relatively fast speed with no direct line of sight. Such mobile channels may 
experience fading due to Doppler shifts caused by the relative speed of vehicles and may be better described 
using Rayleigh fading channel model. 

To decode LDPC encoded codewords transmitted over Rayleigh fading channel, we propose 
modifications to inversion function (5): 


Ek = Xk: Yr — (hg Ees —1)+ a) (7) 


where h,’s are random variables with Rayleigh probability distribution function and are proportional to the 
channel coefficients, described in sub section 2.1. And q,’s are Gaussian random variables with zero mean 
and variance of a*Ny/2, where 0 < a < 1. The perturbation coefficients h, and qx are selected such that 
they exhibit and track the characteristics of a Rayleigh fading channel as well as added receiver noise, 
respectively. However, while it is assumed that channel state information (CSI) is known during the decoding 
process, for practical purposes, h,’s is kept relatively constant across the transmitted frames used in our 
simulations. Channel estimation in modern cellular communication systems (4G and 5G) is typically 
performed at OFDM symbol level or subframe level [22], [23]. This assumption also holds true for V2V 
applications where number of demodulation reference symbols (DMRS) are increased from 3 to 4 to better 
adapt to fast channel variations due to speed of vehicles [24], [25]. Even with the advances in physical layers 
and channel estimation techniques for mobile applications, changing perturbation coefficients at transmitted 
bit level is extremely challenging [26]. Therefore, we found it reasonable to maintain perturbation 
coefficients relatively constant across our simulated frames. For our simulations, the parity check matrix H 
was a 21x21 square matrix corresponding with type 1 projective geometry LDPC (PG-LDPC) code (21,11) 
with row and column weight of 5. Incorporating the negative sign in (7) ensured that the objective function is 
maximized when all parity check equations are satisfied along with the maximized cross correlation of the 
transmitted bit and soft decision received bit. The perturbation coefficients add to the complexity of the 
inversion function as it was originally introduced in (5). Appropriately selecting values of hę and qp is 
extremely important in order to optimize performance of the decoder. In section 3 we will present the results 
of our complexity and sensitivity analysis on how different values of hg and qg impact the performance of 
the proposed decoder in terms of rate of convergence. 

Our proposed decoding algorithm can perform single bit flipping, multi bit flipping and hybrid 
(mode switching) bit flipping. The second modification we made was to the manner by which mode 
switching occurs. Our proposed mode switching method is as following; Start the decoding process by 
initially flipping m bits. As the number of iterations becomes half of the maximum number of iterations Imax; 
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lower the number of bits to be flipped to = And finally, once all individual inversion functions E,,'s become 


positive, switch to single bit flipping (m=1). This method ensures the large steps, associated with multi bit 
flipping, do not miss the desired local maximum, as a known characteristic of multi bit flipping [5], [6], and 
at the same time ensures the rate of convergence is faster than it is of single bit flipping method. The 
proposed decoding process: 


Improved NGDBF algorithm 
Step 1: Compute all syndrome components in s = (S1, S2, =, SJ) = Z. H" (1). If all s;’s are zero output z as 
correctly received codeword, otherwise move on to next step (step 2). 


Step 2: Determine the inversion function Ekg = Xk. Yg — (hk Ye (25> — 1) + qx). Initiate decoding by 
j 


flipping m bits whose corresponding inversion function is the least. 
Step 3: Put the flipped bits in step 2 aside so they do not get flipped in the next immediate iteration. 


Imax 
2 
the individual E's are no longer negative then switch to single bit flipping mode, i.e., flip the bit 


whose inverstion function is the least. Otherwise keep multi bit flipping mode. 
Step 4: Repeat the process until all syndromes are satisfied or the number of iterations has reached Imax- 


Step 4: If the current number of iteration J is equal to then lower the number of bits to be flipped to = If 


3. SIMULATION RESULTS AND DISCUSSION 
3.1. Channel model 

To simulate a wireless communication use case, where codewords are transmitted over Rayleigh 
fading channel, a V2V communication scenario was envisioned where relative speed of approaching vehicles 
varied from 1 mi/h (0.5 m/s) to 120 mi/h (53.6 m/s) corresponding to Doppler frequencies f of 10 Hz to 
1,000 Hz, respectively. Table 1 shows the mean of all Rayleigh distributed random variables used in our 
simulations as channel fading coefficients. In all scenarios, mean of the generated Rayleigh distributed 
coefficients are approximately zero or very close to zero as it should be [18]. Figure 1 also shows the envelop 
of the signal level (transmitted over Rayleigh channel) corresponding to Doppler shifts of 10 Hz, 100 Hz and 
1,000 Hz. As Doppler frequency (shift) increases the level crossing rate also increases which demonstrates 
slow to fast fading characteristics of the Rayleigh channel [27]. For clarity, only one realization of the signal 
envelopes (211 samples) corresponding to Doppler frequencies of 10 Hz, 100 Hz, and 1,000 Hz are shown in 
Figure 1. 


Table 1. Mean of generated Rayleigh distributed channel coefficients 
Doppler Shift (Hz) 10 100 1,000 
Mean of fading coefficients 0.003 0.010 0.034 


a 
2 
i) 
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Figure 1. Signal amplitudes in three scenarios corresponding to Doppler frequencies of 10 Hz, 100 Hz, and 
1,000 Hz 
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Observing the fading rates displayed in Figure 1, we noticed that when the fading rate is the slowest 
(corresponding to Doppler frequency f,, of 10 Hz) the fade values change at such a slow rate that the signal 
level is relatively constant across small enough transmit time intervals. Assuming channel state information 
is known, channel coefficients become highly correlated, and the channel is practically Gaussian [28]. The 
other extreme is when the fading rate is the fastest (corresponding to Doppler frequency fm of 1,000 Hz), and 
the fading is sufficiently fast for the adjacent symbols to be almost uncorrelated. As the Doppler frequency 
increases the amount of correlation (of the received symbols) transitions from highly correlated to almost 
fully uncorrelated. To quantify this phenomenon, normalized Doppler frequency fmT is usually used where T 
is the symbol rate [28]. Obviously, where the transition (from fully correlated signals to uncorrelated ones) 
occurs depends on the system design, the symbol rate, and the depth of interleaving of the code bits. For our 
analysis, we performed the decoding at Doppler frequencies that are decade apart, i.e., 10 Hz, 100 Hz and 
1,000 Hz, respectively. 


3.2. Characteristics and performance of the proposed decoder 
3.2.1. Complexity analysis of the perturbation coefficients h, and qx 

Our investigation has shown that rate of convergence (to the actual solution, i.e., correctly decoded 
codeword) is a function of the perturbation coefficients hę and q,. In other words, decoding performance is 
sensitive to the choice of both hę and qg. To demonstrate this sensitivity, we first kept hg = 1 and varied qxp. 
Figure 2 shows the rate of convergence of the proposed decoder as a function of qg. The case with q, = 0 
(blue trace in Figure 2) is the non-perturbed or legacy case and is considered the baseline. It can be seen that 
the decoder has the fastest rate of convergence when qg exhibits the same characteristics of the added noise. 
This finding is in good agreement with the finding in [6] as well. On the other hand, if q,'s are not properly 
chosen the performance of the decoder is even worse than the legacy (non-perturbed) decoder. The optimum 
dx did not show any sensitivity to Rayleigh channel coefficients (or Doppler frequency) nor to the signal to 
noise ratio. Figure 2 shows the rate of convergence for the worst case signal to noise ratio in our study that is 
SNR=0 dB. 
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Figure 2. Rate of convergence as a function of perturbation coefficient q,, SNR = 0 dB 


The next step was to determine the value for the second perturbation coefficient h, that optimizes 
the decoder performance. To do so, we kept the coefficient q, the random variable drawn from Gaussian 
distribution that produced the fastest rate of convergence shown in Figure 2 and then varied the hę in 
multiple scenarios with different Doppler frequencies. Figures 3 and 4 show the rate of convergence as a 
function of perturbation coefficient hg where Doppler frequency was 10 Hz, and 1,000 Hz, respectively with 
no measurable difference between Doppler frequencies of 10 Hz and 100 Hz. It can be seen that the decoder 
has the fastest rate of convergence when hę is approximately equal to the channel coefficient Rẹ in both 
scenarios. The perturbation coefficients h, did not demonstrate any measurable sensitivity to SNR. The 
proposed decoder’s performance is shown for the worst case signal to noise ratio in our study that is 
SNR=0 dB. 


Improved noisy gradient descent bit-flipping algorithm over Rayleigh fading channel (Reza Biazaran) 


2706 O ISSN: 2088-8708 


Our introduced perturbation coefficients q, and h,, when properly selected, improve the 
performance of the decoder as compared to the non-perturbed or legacy decoder and as shown in Figures 2 
through 4. However, the improvement comes at a cost and that is both the receiver noise variance and 
channel coefficients have to be known by estimators that are external to the decoder, hence added complexity 
to the decoder and receiver design. With the added complexity in mind, our proposed decoder is better suited 
for applications such as V2V, described in section 2, where capability to estimate noise variance and channel 
coefficients is assumed. 


Objective Function 


0 1 2 3 4 5 6 7 8 9 101112131415 


Number of iterations 


Figure 3. Rate of convergence as a function of perturbation coefficient hg, SNR = 0 dB, fm = 10 Hz 
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Figure 4. Rate of convergence as a function of perturbation coefficient hg, SNR = 0 dB, f = 1,000 Hz 


3.2.2. Error performance of the proposed decoder 

Figure 5 shows the error performance of our proposed decoder comparing to NGDBF decoder while 
codewords are transmitted over Rayleigh fading channel with Doppler frequency of 10 Hz. Our proposed 
decoder outperforms the legacy NGDBF algorithm across the full range of measured signal to noise ratios. 
Another performance metric, that becomes especially important in use cases with low latency requirements, 
is the decoder’s rate of convergence since it directly relates to the overall processing delay and ultimately end 
to end latency. Figure 6 depicts the average number of iterations over hundred frames to achieve the error 
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performance of our proposed decoder shown in Figure 5. Again, our proposed decoder outperforms the 
legacy NGDBF decoder across the full range of measured signal to noise ratios. On average, our proposed 
decoder shows 19.98% improvement in average number of iterations required to obtain the said error 
performance. 
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Figure 5. BER curves of the improved and legacy NGDBF decoders over 100 frames, 
max iterations: 15, fm = 10 Hz 


14 
12 
10 NGDBF 
(AWGN) 
8 Decoder 


— - + This Decoder 


Average Number of Iterations 
(ep) 


4 
2 
0 
0 1 2 3 4 5 6 
Eb/NO (dB) 


Figure 6. Average rate of convergence for the improved and legacy NGDBF decoders over 100 frames, 
Max iteration: 15, fn = 10 Hz 


There was no measurable difference in the performance of the decoders over Rayleigh fading 
channel with fm = 100 Hz. Our decoder’s performance, however, was improved in case of fm = 1,000 Hz 
as it can be seen in Figure 7. The improvement should be expected since LDPC codes are random error 
correcting codes, and they perform best if introduced errors into the codewords are random in nature and 
statistically independent. As the Doppler frequency increases the fading rate becomes sufficiently fast so the 
adjacent symbols are uncorrelated, and the errors appear random. Our finding agrees with the studies in [11] 
and [28] as well. Figure 8 depicts the average number of iterations over hundred frames to obtain the error 
performance of our proposed decoder shown in Figure 7. Again, our proposed decoder outperforms the 
legacy NGDBF decoder across the full range of measured signal to noise ratios. On average, our proposed 
decoder shows 16.87% improvement in average number of iterations required to obtain the said error 
performance. 
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Figure 7. BER curves of the improved and legacy NGDBF decoders over 100 frames, max iterations: 15, 
fm = 1,000 Hz 
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Figure 8. Average rate of convergence for the improved and legacy NGDBF decoders over 100 frames, max 
iteration: 15, fn = 1,000 Hz 


The results of our simulation, can also be analyzed from an information theory perspective. 
Codewords that are transmitted over a Rayleigh fading channel may experience block error (in case of slow 
fading) when transmitted frames experience deep fades, or they can undergo a fast fading when introduced 
errors are less correlated and error patterns are more random in nature. In case of slow fading, corresponding 
with Doppler frequencies of 10 Hz and 100 Hz in our simulations, the outage probability can be 
approximated when knowledge of channel state information is available at the receiver. The outage 
probability in this case has been shown to be inversely proportional to signal to noise ratio at a given 
transmission rate [29]—[31]. It is equivalent to the fact that the effective channel capacity increases as signal 
to noise ration increases. Our results shown in Figures 5 and 7 also show the same trend and agree well with 
established results in [29] and [30]. In case of fast fading scenario, corresponding with Doppler frequency of 
1,000 Hz in our simulation, channel capacity has been shown to be directly proportional to signal to noise 
ratio [29] as well and the error probability decreases as the signal to noise ratio increases. Additionally, [28] 
has shown that the error performance improves as channel transitions from slow fading to fast fading. Our 
proposed decoder’s performance improves by a factor of 10 (probability of bit error of 1x107” at SNR of 
5.5 dB in Figure 7) when fast fading is experienced comparing to probability of bit error of 1x10~* at the 
same SNR as shown in Figure 5 when slow fading scenario is experienced. Again, our results are in very 
good agreement with results in [28]. 
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4. CONCLUSION 

In this paper, we proposed an improved noisy gradient descent bit flipping decoder that performs 
well over Rayleigh fading channel that is a more realistic channel model for a lot of modern applications with 
low latency requirement. We demonstrated that our proposed decoder outperformed both non-perturbed and 
legacy NGDBF decoders that are designed for AWGN channels. Measurable improvements are observed in 
terms of error performance and processing delay of the decoder. Additionally, our proposed decoder 
outperformed legacy NGDBF decoder across a wide range of mobile channel conditions quantified in terms 
of Doppler frequency to simulate slow and fast fading channels. The fact that our decoder’s performance 
improved as fade rate increased agrees very well with established results reported in other studies. Our 
empirical results also prove the hypothesis that the decoder performance is optimum (or near optimum) when 
the perturbation coefficients exhibit and track the statistical characteristics of the channel coefficients. 
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